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BROWNIAN COAGULATION AND A VERSION OF 
SMOLUCHOWSKI'S EQUATION ON THE CIRCLE 

By Ines Armendariz 

Universidad de San Andres 

We introduce a one-dimensional stochastic system where parti- 
cles perform independent diffusions and interact through pairwise 
coagulation events, which occur at a nontrivial rate upon collision. 
Under appropriate conditions on the diffusion coefficients, the co- 
agulation rates and the initial distribution of particles, we derive a 
spatially inhomogeneous version of the mass flow equation as the par- 
ticle number tends to infinity. The mass flow equation is in one-to-one 
correspondence with Smoluchowski's coagulation equation. We prove 
uniqueness for this equation in a broad class of solutions, to which 
the weak limit of the stochastic system is shown to belong. 

1. Introduction. Coagulation models describe the dynamics of cluster 
growth. Particles carrying different masses move freely through space, and 
every time any two of them get sufficiently close there is some chance that 
they coagulate into a single particle, which will be charged with the sum of 
the masses of the original pair. 

In 1916, Smoluchowski [18] considered the model of Brownian particles 
moving independently in three-dimensional space, such that any pair coag- 
ulates into one particle upon collision. He derived a system of equations, 
known as Smoluchowski's coagulation equations, that describes the time 
evolution of the average concentration /if (to) of particles carrying a given 
mass m = 1,2,.... In this original work, Smoluchowski ignored the effect 
of spatial fluctuations in the mass concentrations, the equations we write 
below are thus a natural extension allowing diffusion in the space variables: 

fi t (m) = ^a(m)A x fi t (x, m) + - ^ re(m', m")[i t (x, m')fi t (x, to") 

m! +m"=m 



Received May 2006; revised July 2009. 

AMS 2000 subject classifications. Primary 60K35; secondary 82C2f . 
Key words and phrases. Coagulating particle systems, hydrodynamic limit, Smolu- 
chowski's equations. 

This is an electronic reprint of the original article published by the 
Institute of Mathematical Statistics in The Annals of Applied Probability, 
2010, Vol. 20, No. 2, 660-695. This reprint differs from the original in pagination 
and typographic detail. 



1 



2 



I. ARMENDARIZ 



— Ht(x, m) k(tti, m')fjLt(x, m'). 
m' 

The differentiated term on the right describes the free motion of a particle 
with attached mass m as a Brownian motion with diffusivity rate a(m). The 
kernel k is determined by physical considerations, it regulates the intensity 
of the coagulation dynamics. The first sum corresponds to the increase in the 
concentration resulting from the coagulation of two particles whose masses 
add up to m. The second sum reflects the decrease caused by the coalescence 
of a particle carrying mass m with any other particle in the system. Coag- 
ulation phenomena have been studied in many fields of the applied science; 
we refer to Aldous's review [1] for a comprehensive survey of the literature. 

This paper is concerned with the approximation of Smoluchowski's equa- 
tions by stochastic particle models. Concretely, we are interested in iden- 
tifying the solution to the coagulation equations as the mass density of a 
system of interacting particles, when the particle number tends to infinity. 
This problem has been much studied in the spatially homogeneous case, both 
for discrete and continuous mass distributions, and with different choices of 
coagulation kernel k (cf. [4, 5, 8, 15] and references therein). The relevant 
stochastic process for these models is the Marcus and Lushnikov process 
[13, 14]; this is the pure jump Markov process where clusters of size m and 
m' coagulate into a single cluster of size m + m' at rate K(m,m'). 

In the spatially inhomogeneous case, on the other hand, the coagulation 
mechanism is highly dependent on the relative position of the particles, hence 
the space dynamics plays a predominant role in the particle interactions of 
the stochastic system. In the original problem proposed by Smoluchowski, 
for instance, pairwise collisions and the ensuing coagulation events are com- 
pletely determined by the Brownian paths. The first result for this model 
was obtained in 1980 by Lang and Xanh [12] for the case of discrete mass 
and constant a, k, in the limit of constant mean free time. No progress was 
made until the forthcoming paper [16], where Norris proves convergence for 
both discrete and continuous mass distributions, and variable coefficients in 
a class that includes the Brownian case. 

Over the past few years there has been considerable interest in spatial 
models with stochastic dynamics of coagulation. In these models, particles 
coagulate at some rate while they remain at less than a prescribed distance. 
Deaconu and Fournier [3] consider the case when this distance is independent 
of the particle number, and let it go to zero after taking the weak limit. The 
moderate limit, where the range of interaction is long in the microscopic 
scale, is studied by Grosskinsky, Klingenberg and Oelschlager in [7] in a 
regime where the dominating particle interaction are shattering collisions. 
In the articles [9, 10], Hammond and Rezakhanlou work in the constant 
mean free time limit, for dimensions d>2. 
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In this paper we introduce a diffusion model where coagulation occurs on 
collision as a result of a random event: N mass-charged particles perform 
independent one-dimensional diffusions, and whenever two particles are at 
the same location they may coagulate at a positive rate in their intersection 
local time. The new particle is assigned the sum of the masses of the incoming 
particles, and the process continues. 

This model is motivated by the problem of establishing the large-scale 
dynamics of a system of Brownian particles confined to a thin tube, inter- 
acting through pairwise coagulation when any two of them get close enough. 
It would be interesting to determine whether, under proper scaling of the 
tube and particle radii in terms of the system size N, the higher-dimensional 
model can be replaced by the simpler one-dimensional one. 

Note that in one dimension the problem of instantaneous coagulation 
on collision is not interesting: due to the recurrence properties of the one- 
dimensional Brownian path, in the limit we would instantaneously see the 
distribution of the total mass in the system among clusters of macroscopic 
size. In fact, in order to keep the average time T during which a tagged par- 
ticle does not undergo a collision constant, it is necessary to set coagulation 
rates which are inversely proportional to the particle number N. The model 
under consideration is thus of the constant mean free time type, and in this 
sense it is a one-dimensional version of the models studied in [9, 10, 12] and 
[16]. 

We treat the case where the mass dependent diffusivity rates blow up as 
the mass goes to zero, combined with our choice of rates, this leads to a 
large-scale model favoring coagulation of large and small particles. 

We describe the particle system and state the main results of the paper 
in Section 2. Section 3 contains a compactness result. The next step in the 
analysis is to prove convergence to a hydrodynamic limit. This is shown to 
verify a spatial version of the mass flow equation, which is closely related 
to Smoluchowski's coagulation equation. This is the content of Section 4. In 
Section 5 we derive a uniqueness result for the solutions to a broad family of 
such equations, thereby obtaining a law of large numbers for the empirical 
processes of the microscopic model. 

2. Notation and results. Consider a positive integer N. Let T stand for 
the one-dimensional torus, and 1R+ for the half line [0, oo). Let 

pN = p* (dx l ,dm 1 ;...;dx N , dm N ) 

be a sequence of measures on (T x ]R + ) Ar which are symmetric on the pairs 
(x l ,m l ) and supported on X^ m * = 1- Denote by M\(T x K + ) the space of 
probability measures on T x R + endowed with the weak topology, and by 
Aif(T x M + ) the space of positive and finite measures on T x R_|_. 
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For a complete separable metric space T, we will denote by D(M + ,T) 
(or D(I,T), I an interval of the real numbers), the set of right continuous 
functions with left limits taking values in T, endowed with the Skorokod 
topology. 



2.1. The particle model. Let 

$(m, m!) :I + xK + — >R + 

be a nonnegative, symmetric kernel with the property that it vanishes when 
either of the coordinates equals 0. 

Given a point {(x*,m 4 )} in (T x R + )-^, define 

(X t To ,M^ ) = (x i ,m l ), l<i<N, 

and set To = 0. Let fcgNU {0}, and suppose that the process 

e = {(JCM)h<i<N G C([0,T k ],T N ) x D([0,T k ],M") 

has already been defined up to the time Tj., a stopping time with respect 
to the <T-algebra J~t = <r{£^,0 < s < t} generated by S, N . Consider then a 
family of ./V independent Brownian motions '*} on T with corresponding 
diffusion coefficients a{NM^ ) and initial positions f3 Q ' 1 = . For each 
pair i < j, denote by L k M the intersection local time of the ith and jth 
particles; that is, the local time at the origin of the difference /3 fe '* — f3 k 'K Let 
{e k, ^}i<i<j<N be a sequence of (^) independent, parameter one-exponential 
random variables, and define the stopping times 

f <|)(.V.\/' . .V.\/.' ) ... , A 
T k+l = mm{^ +1 }, Til, = infji > 0, ^ %L T -^L^ > ^ j. 

Let then Tfe +1 = + T^+i, and for 1 < i < N, set 
X\ = P t -T h i T k <t< T k+ i , 
M\ = M^ T k <t<T k+1 , 



M l Th + , if f k+1 = T l k 3 +V for some j > i; 
0, if T fc+ i = Tl % +1 , for some j < i; 



ML , otherwise. 



We will denote by U 3 the intersection local time of the ith and jth parti- 
cles X 1 and X J . Note that between two consecutive stopping times T k and 
T k+1 the identity L ij = L k ^ holds. 

The dynamics is well defined except for those configurations where two or 
more coagulation events occur simultaneously. Let us briefly show that the 
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set of such configurations has measure zero and can therefore be neglected. 
We first show that this is the case on the time interval [0, Ti], when X 1 = 
1 < i < N. The ith. and jth. masses will coagulate at a time belonging 
to the support of the measure dL l i , which equals the zero set of /3 1 '* — /3 lj . 
Fix two pairs of indices i < j and k <l. If the four indices are different, then 
0i,i _ 0i,j anc j 0i,k _ p er form independent diffusions, and from the fact 
that point sets are polar for Brownian motion in 2 or higher dimensions it 
follows that these diffusions do no vanish at the same time. Let there be a 
repeated index: j = k, say. Set a = ~^~^r^fjy~ ' ■ Then 

U t = a \f}¥ ft* and V t = fti-Pl> 1 

are independent diffusions in T 2 . The previous argument implies that with 
probability 1 they never vanish simultaneously, and hence the same applies 
to /3 1 '* — /3 lj and /3 lj — We conclude that in any case the zero sets of 
0i,i _ anc j pi,k _ are ciisj oint , with probability 1. By repeating the 
argument on each interval [Tfc,Tfe+i],fc > 1, it follows that outside a set of 
measure zero there are no conflicting coagulation events. 

2.2. Martingales. Let P N be the measure on C(R+,T N ) x £>(R + ,R^) 
determined by the process £ N . There is the representation 

M l t =m i + f M j dE ij - j dE ki , 

i<j<N l<k<i 

where dE %] is a counting measure, £V([0,t]) = 0,l with P^-probability 1, 
depending on whether the ith and jth particles have coagulated by time 
t > 0, and 

H?=E<HM)- [' mM - N NMi) dL* 

J 

is a martingale. 

In general, if / is a bounded function on (T x with two continuous, 

bounded derivatives in the space coordinates, then 

m)- / x>(jva£)— k(&)ds 

Jo ^ dX^ 

-/'E[/(g)-/fc)] "" Ml> ^ 

is an {Tt, P)-martingale. Given £ = {(X 1 , M 1 )} , here is given by 

(J^AP + AP), if k = i; 
(C 3 ) k ={ (X3,0), if k = j; 

(X k ,M k ), otherwise. 
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2.3. Scaling. The reason for the choice of scaling of the coagulation rate 
is quite straighforward: if a hydrodynamic description is to hold, then it is 
necessary that O(N) mass charged particles remain in the system at all times 
(note that although total mass is conserved, the number of particles carry- 
ing positive mass decreases by one after each coagulation event). Therefore 
a generic particle will see some fraction of the other O(N) mass charged 
particles over any fixed time interval, while it would still be expected to 
coagulate with only 0(1) of them. This forces the rate $ to be typically of 
order 1/iV. 

2.4. Assumptions. We will consider the mapping 

U N :C(R+,T N ) x D(R+,M^) -+D(R+,Mi(T x R + )) 

such that 

U N ({X i ,M i })= ]T MH^i^Mi) 

\<i<N 

and denote by Q N the measure on D(M + , A^i(T x M+)) induced by ITv, 

q" = p" on- 1 . 

In order to derive a hydrodynamic limit for Q , we need to specify some 
technical conditions on the coalescing kernel the diffusion coefficients a, 
the initial measure and the profile v. 

The kernel satisfies a Lipschitz condition away from the origin: for each 
L > there exists a positive constant T(L) such that 

(2.1a) |$(m + to,", to') — $(m, m')\ < T(L)m" whenever to > L. 

It will also be assumed that there exists < p < 1/2 such that 

(2.1b) $(m,m') < c(m p + m' p )l [m>0 . m , >0] 

for some positive constant c. 

We assume that the mapping a : (0, oo) —> (0, oo) is nonincreasing, so that 
particles diffuse at a slower rate as they gain mass. As a consequence the 
kernel $ and the diffusion coefficients do not grow simultaneously. We set 
a(0) = 0. 

We are ready to introduce the coagulation rates k : M+ — > R.+ appearing 
in the hydrodynamic equation, 

k(to, to') = $(to, m')[a(m) + a(m )]. 

In order to study convergence and derive the uniqueness of the limit it will 
be useful to consider 



uj(m) = [1 + c + a(l)] [to p + a(m) + 1], 
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it verifies K(m,m') < ui(m)uj(m'). We will then require that 
(2.2) a(m)~ 1 / 2 o;(m) be a subadditive function of m. 

Conditions (2.1a), (2.1b) and (2.2) are for instance satisfied by 

<]?(m, m) = C(m a + m a ) and a(m) = — ^l[ m >o] with a < — and j3 < 1. 



77? 



In this case K(m, m) = C(m a + m a )(^- + ^g- )• 



We will assume that there exists an initial profile v € A4\(T x M + ) such 
that the empirical distributions ^2 i m l 5( x i Nm i} converge in distribution to 
by as N — > oo, where by St x i tNm i) (resp., 5 U ) we denote the Dirac measure 
with a unit atom at (x l ,Nm l ) (resp., v). 

The initial measures P$ will satisfy 



(2.3) E p o N [nJ2(^? 



<C and E F ° 



for some constants C, C > 0, uniformly in N . In particular both {m, v) and 
(a(m) 2 /m,u) are finite. In fact, the following assumption will hold: there 
exists a finite measure v*{dm) such that 

/ 0(771)^ ^ 

(2.4) v{dx,dm) <v*(dm)dx with ( m -\ ,/>*>< 00. 

\ 771 



2.5. Results. The first theorem of the paper is a tightness result: 



Theorem 1. Assume that (2.1a), (2.1b), (2.2) and (2.3) hold. Then the 
sequence of measures Q N on D(R + ,A4\(T x is relatively compact, and 
all limit points are concentrated on continuous paths. 



The next two results concern the properties satisfied by any weak limit 
of the empirical distributions as we pass to the limit in the particle number. 
The first result provides some estimates that will ensure the hydrodynamic 
equation is well defined, then Theorem 2 identifies this equation, thereby 
establishing an existence result. 

Given a kernel admitting a representation fj,(x, dm) dx S A4f(T x ]R + ) and 
a bounded test function f(m), we will denote by ((/, fi)) the single integral 
f (m) n(x , dm) . This clearly determines a signed measure on T by 



I h(x)((f, fi}} dx = / h(x)f(m)n(x,dm)dx. 

T JTxR + 
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Proposition 1. Let Q be a weak limit of the sequence Q N of probability 
measures on C(M> + ,Mi(T x Then Q is supported on the set of paths 

fj,.(dx,dm) whose marginal fit(dx,M + ) <C dx on T for all t, /j,t(dx,dm) = 
Vt(x,dm)dx. Moreover, the following inequalities hold with Q -probability 1: 



(2.5) 
and 



sup 

t>o 



;(m) 



m 



, v f 



< oo 



(2.6) 

for any fixed final time T . 



sup(m, fit) < oo 

t<T 



Denote by C fe 2 (T x R+) the space of continuous, bounded functions on 
T x R + which have continuous, bounded derivatives in the space variable up 
to the second order. 

Theorem 2. Let Q be a weak limit of the sequence Q N , as in Proposi- 
tion 1, and consider f in C^(T x R+). Then, with Q -probability 1, a path 
fj> s (dx,dm) satisfies 

/■*/! d 2 f \ 



(2.7) 



+ 



[f(x,m + m') - f(x,m)] 



yT7M + xR + 771 

x v s (x, dm)v s (x, dm') dx ds, 



n{m, m) 



if t>0. In this equation v is the initial profile of the model, and the kernel 
v s (x,dm) is such that 

H s (dx,dm) = v s (x,dm) dx for alls >0, Q-a.e. 

Equation (2.7) describes the evolution in time of the mass flow: if we de- 
compose its solutions as v s (x,dm) dx = mv s (x,dm) dx, an elementary com- 
putation proves that v dx satisfies Smoluchowski's coagulation equation with 
kernel k. Theorem 1 asserts that all weak limits of the measures Q N are sup- 
ported on A^i(T x M+); in terms of the concentration densities v s (x,dm), 
this means that mass is conserved, or equivalently, that there is no gelation 
phenomenon. 

The method applied to derive these results relies heavily on stochastic 
calculus computations, we try to make these quite detailed in the proof of 
Theorem 1 and give only an outline later on. 
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In [16], Norris introduces a method for proving existence and uniqueness 
for a general class of <i-dimensional diffusion models with coagulation; briefly 
put, this consists on approximating the corresponding version of (2.7) in his 
paper by a system that depends on the coalescing kernel k only through its 
values on a given compact set. In Section 5 we develop a simplified version 
of his technique to obtain a uniqueness result for the solutions of a broad 
family of mass flow equations. Section 5 may be read independently of the 
rest of the paper. 

Some brief consideration shows that the right-hand side of (2.7) is well 
defined in a proper subset of C(R+, Mi(T x R+)) consisting of those paths rj 
whose marginal rjt(dx,M + ) has a density with respect to Lebesgue measure 
satisfying some integrability conditions. It is easy to see that in fact (2.5) 
and (2.6) are enough, and then Proposition 1 says that all weak limits of the 
sequence Q N are supported on configurations where the right-hand side of 
(2.7) can be evaluated. This observation motivates the following definition: 
we will denote by V(lo) the subset of C(R+, «M/(T x R+)) of those paths r\ 
whose marginal r]t(dx,M + ) <C dx, and such that (2.5) and (2.6) are satisfied. 
Note that the map oj depends on the diffusivity a and the coagulation rate 
<]?, and then so does T>(uS). 

As a particular case of Theorem 3 in Section 5, we have: 

Corollary 1. Assume conditions (2.1b), (2.2) and (2.4) on the coag- 
ulation rate k and the initial measure v , respectively. Then for any T > 0, 
(2.7) has at most one solution {fJ>t}o<t<T in T>(u). 

The four preceding results imply that the sequence of probability measures 
Q N converges to the Dirac measure concentrated on the unique solution in 
V{u) to (2.7). 

3. Existence of a weak limit. In order to simplify notation, we will often 
omit the dependence of the masses and positions on the time parameters 
whenever we think that this would not lead to confusion. For instance, in 
an integral where time is parametrized by s, M % and X % should be read as 
Mj and X l s , respectively. 

Throughout the article, T will denote a positive constant. Unless we are 
particularly interested in keeping track of its growth or dependence on the 
parameters, we will use the same letter V to denote constants on consecutive 
lines which may be different, or constants appearing in totally unrelated 
computations. 

Let us consider a fixed final time T > for the rest of the paper. We will 
prove a version of Theorems 1, 2 and Proposition 1 on the compact interval 
[0,T]; the fact that the value of T is arbitrary will then imply that these 
results hold as stated in the previous section. 
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The following estimates will be necessary to derive Theorem 1; we post- 
pone their proofs until the end of this section. 

Lemma 1. There exist nonnegative constants C(T),C'(T) which depend 
on the diffusivity a, the kernel 4> and the bounds appearing in (2.3), such 
that 



(3.1) 



(3.2) 



E 1 



E pN NY,[M l T } 2 

i 

M i M j ^{NM\NM j ) dL ij 



i<j 



< C(T) and 

< C'{T) 



hold uniformly in N . 

Lemma 2. Given e > 0, there is 5 > such that 



lim P N 



sup / V M i M J '\>(.\ M : . NM j ) dL ij > e 
s<t<rJs 



o< 

t-s<S 



< e. 



Denote by C&(T x M + ) the space of bounded, continuous functions on 
T x IR + with the topology determined by uniform convergence over compact 
sets. Let {fk,k G N} be a dense, countable family in C&(T x R+). Then the 
distance 

, 1 \{fk,»)-{fk,v)\ 

^-^l + \(f k ,^)-(f k ,u)\ 

defines a metric on A4i(Y x M + ) which is compatible with the weak topology. 
There is the associated modulus of continuity 

cj m (7) = sup Q(n t ,(j, s ). 

0<s<t<T 
t— s<7 

Proof of Theorem 1. We refer to Chapter 4 in [11] for a presentation 
of the Skorokod's topology as well as the characterization of the relatively 
compact sets in D([0, T], A^i(T x M+)). Note that condition (ii) below im- 
plies that, provided the sequence Q N has limit points, these will be sup- 
ported on C([0,T],Mi(T x R+)). 

By a version of Prokhorov's theorem applied to this setting (cf. [2], Chap- 
ter 3), the theorem will follow if we can show that: 



(i) 



For every e > 



lim lim sup Q 



N 



sup fjbt(m > M) > e 

0<t<T 
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and 
(ii) 



For every e > limlimsupQ [^(7) > s] = 0. 

74.0 JV-VOO 



Note that (m,^) is nondecreasing for < t < T, Q fact we will 

repeatedly use in the course of the article. Then (i) is an easy consequence 
of (3.1) in Lemma 1 and Chebyshev's inequality. 

In order to conclude (ii) it will be enough to prove that given / G Cj (T x 
/ Lipschitz in m, we can control 



Q 



sup I (/,//*) - (/,Ms)| >e 

0<s<t<T 
t— s<7 



< £ 



provided iV and 7 are taken to be sufficiently large and small, respectively. 
Define the stopping time 



t = m£\ t > 0,max(M|) 



1 + y^m j a(Nm j ) 



> 



JV-l/4 



by Chebyshev's inequality we compute 



[c(r)(i + c")]V 2 

ATl/4 



where C(T) and C are the constants appearing on the right of (3.1) and 
the second inequality in (2.3), respectively. By stopping the process as soon 
as t is achieved, we may assume that 



(3.3) 



P 



N 



max(iW| 



1 + m? a{Nm^ 



< AT" 1 / 4 



Applying Ito's formula to /, we have 
= f^£ j M i ^-{X\NM i )dX i 

Js „■ OX 



+ \ I J2 Mi p2( xi > NM >( NMi ) ds 

Js i x 

+ f ^2[F N (X i ,M i + M j )- F N {X i ,M i )- F N {X\M j ))dE i3 , 

J s 



where we denote F^{x,m) = mf(x,Nm). 
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Let 7 > 0. Doob's inequality, (3.3) and the monotonicity of a imply 



P 



N 



sup 

0<s<t<T 
t— s<7 



>iV 



sup 

0<f<T 



< P 



[ ^M i ^(X\NM i )dX i 

Js i x 

f V M i {X\NhP) dX i 
Jo Ox 



e 

> 6 



N 1 / 4 ' 

where T is a positive constant that does not depend on N. By taking 7 such 
that 

(3.4) 

we also obtain 



c 



d 2 f 



dx 2 



7 <f! 

, 3 



sup 

.0<s<t<T 
t— s<7 



J S ^ (X i ,NM i )a(NM i ) ds 



> 



< 



It remains to estimate the Poisson integral 



I Y^i F N{X i ,M i + M j )-F N {X i ,M i )-F N (X\M j )]dE lj 

J s 



=n F (o,t) -n F (o,s) + i F (s,t), 

if I F (s,t) denotes the integral 

/ Yl ^ Fn ' Mi + Mi ) ~ Fjv ( xi ' Mi ) ~ Fjv ' Mi ■ 



$ (JVM*, .WAP) 
' iV 



and 7i F (0,t) is a martingale collecting the remaining terms. Its quadratic 
variation is given by 

*t 



[ J2 [Fn {X i , AT + M j )-F N {X\M i )- F N {X i ,M j )} 2 



$(NM\NM3) 
N 



dVK 



i<3 

Note that 

\Fj^(x,m + m ) — Fj^(x,m) — Fjy(x,m )| < r(/)[(m + m ) A (JVrara')]. 

In particular, due to assumption (3.3) on the mass sizes, (3.2) and Doob's 
inequality, we obtain 



P 



N 



sup \n F (s,t)\>- 

0<s<t<T D 
t— s<7 



< 



r(/,e,r) 

ATl/4 ' 
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which will decay to as we pass to the limit N — > oo . Finally, 

|!>(«,*)|<r(/) f ^M i M j <S>{NM\NM j )dL ij . 
Js i<3 

The result now follows by taking 

7 < 7i A 5, 

where 71 satisfies (3.4) and 5 is the value given by Lemma 2 when e is set 
equal to e/6[l + □ 

Proof of Lemma 1. The proof of this lemma will follow from repeated 
applications of Ito-Tanaka's theorem; see [6] for an exposition of this and 
related formulas. We write 

(3.5) iVV[M t i ] 2 = AfV[m i ] 2 + H i + 2 I VM]M^(iV 'M l s ,N 'Alj) dL i3 , 

where Tit is the P^-martingale 

U t = f V 2NMM \dE- - mS*™*) d L« 
Jot?, L N 

We focus on the last term of (3.5). Given ( > 0, let g c E C(T) n C 2 (T - {0}) 
be a positive, even function that equals |x| in a small interval containing the 
origin, vanishes outside [—1/4,1/4] and satisfies sup x£ j g^(x) < (. We then 
have 

/ M*M j <I> (NAP , NAP ) dl. IJ 

(3.6) 

= A l (0, t) - A 2 (0, t) - A 3 (0, t) - A 4 (0, t) - As(0, t) 



with 



A 1 (0,t)= MiMi$(NMi,NMi)g c (Xi-X() 

l<i<j<N 

- rn i mP^(Nm\Nm j )g (: (X i - X j ), 

A 2 {0,t)= f Y M i M j ^(NM i ,NM j )g'^(X i — X j )[dX i — dX j ] 

l<i<j<N 
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i r l 

A 3 (0,t) = - / S"M l M^{NM\NM 3 )g'l{X l - X^) 

x [a(NM i ) + a{NM j )\du. 

The function g" appearing in the formula for ^3(0, t) stands for what is left of 
the second derivative of (in the sense of distributions) after substracting 
25q, 5o the Dirac measure at the origin. The terms A^(0,t) and A^(0,t) 
correspond to the coagulation martingale and its compensator, 



A 4 (0,t) = [ t J2D N (M i ,M*,M k )g c (X i - X*) 
Jo 7~,. 



i<k 
3 



dE 



lk ${NM\NM k ) dl}k 



N 



M0,t)= fY. D N{M\Mi,M k )g c {X*-Xi f {NMl N NMk) dV\ 
•*° i<k 
5 

where D^(m,m' ,m") is defined as 

(m + m")m'<&(Nm + Nm" , Nm') 

— mm'$>(Nm, Nm') — m"m'&(Nm" , Nm'). 

In deriving the formula for A4, we have used that at the time when the 
masses M l and M k coagulate, the ith and feth particles are occupying the 
same position. 

We will study these terms separately. Replacing <J>(m,m') < c[m p + m /p ] 
in the definition of A\ gives 

(3.7) |A x (0,t)| <4cCX)Afj[iVAfj]*'. 

i 

The bounded variation term A^(0,t) may be similarly controlled, 



(3.8) 



|A3(0,t)|<r(C) y l + Y,Ml[NMif l + £Wa(W) 



ds. 



In order to bound A5, we first notice that by the Lipschitz assumption 
(2.1a) on <I> we have 



(3.9) 



\<£>(Nm + Nm",Nm) -$(Nm,Nm')\ 

< T(l)Nm"l {Nm > 1} + (1 + [Nm'}* + [Wf)l {JVm<1} . 
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\A 5 (0,t)\<T( 



(3.10) 



V M i M k &(NM i , NM k ) dL ik 
i<k 

+ (l + J2Mi[NMi] 
j 

ft 1 

x / ^—^(NM\NM k )dL 

J ■ , u 



i<k 



We have 



rt i 

/ Y T72®(N M\NM k )dL tk 

i 

T.- 

^ N 



N' 
t 



i<k 



dE *_WM>,NM*) iLit 



The last term above is a martingale, hence 



E 1 



I S"\®( NM \ NMk ) dLik 

J ° i<k 



< 1 



and 



E 1 



< 



J Yljf2®( NMi i NMk } dIjlk 

ft 1 

/ J2^3^ NMi ' NM ^ dLik 



1 + B 



< 2. 



We take expectations in (3.6) and combine with (3.7), (3.8) and (3.10) to 
obtain 



(1-TQE 



< 4c(E } 



^M i M k <S>(NM\NM k ) dL ik 

eft 



i<k 



+ T(C) J* E pN (l + J^M^NM^A (l + J^m^iNm 1 ) 



ds 
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i ' ^° i<k 

i 

1 f* 

l + ^m^lW) 2 J E pN 1 + 



ds 



I? 



In order to derive this last bound we have used Holder's inequality and the 
fact that p < 1/2. From (2.3), we conclude that 



(3.11) 



(i-rc)£ p 
<r'c 



l + E 



M i M k <S>(NM\NM k ) dV k 



i<k 
pN 



rt r 

+ r(c) / E pN l + ivViM^ 2 

JO : 



ds. 



Choose ( < l/(4[r + r']), where T and T' are the constants appearing in 
the first line and in front of the first term on the right above, respectively. 
Combining (3.11) with (3.5) we get 



<r 



E 1 



1 + N^2[m 1 } 2 + J (l + Nj^lM*} 2 ) ds 



Estimate (3.1) now follows from Gronwall's lemma and conditions (2.3) on 
the initial distribution of masses, and (3.2) is immediate from (3.11). □ 

PROOF of Lemma 2. Choose ( > and 5 > such that 

2 2 

4C[l + r][l + C + C(T)]<g and 4*T(C)[1 + C(T)][1 + C] < ^, 

where T and T(Q are the constants appearing on the right of (3.10) and 
(3.8), respectively, C and C are the constants from assumption (2.3), and 
C(T) the bound established in Lemma 1. Set the parameter of g equal to a 
value of C determined as above. 

As in the proof of Theorem 1, we will stop the process at the finite stopping 
time r A A T, where r is the stopping time defined in the proof of Theorem 
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IT 



1, and 



r c = inf|t>0,iV^[M^ 2 >i|. 



By Lemma 1 and the choice of £, we have 

P N h<T]<C(T)C<^ 

if e is small enough. We will thus assume that P N is supported on 

(3.12) max{Afj} 



The proof will now follow by estimating the variation of the terms {^4i}i<i<5 
on the right of (3.6). Let Ai(s,t) = Ai(0,t) - Ai(0,s),l< i < 5. 

By Chebyshev's inequality, (3.7), (3.8), and the choice of £, 5, we have 



(3.13) 

(3.14) 

and 
(3.15) 



P 



N 



P 



N 



sup |Ai(s,t)| > - 

,0<s<t<T 0. 
t-s<6 

sup \A 3 (s,t)\>^ 

,0<s<t<T 0. 
t-s<5 



sup \A 5 (s,t)\>- 

0<s<t<T 0. 
t-s<<5 



< 



10' 



e 

< — 
~ 10 



e 

< — . 
~ 10 



The quadratic variations Q2(0,t) and Q^(0,t) of the martingale terms A2 
and A4 satisfy 



Q 2 (0,t)<T(C,T) max{M^} 



< r(c,T) 

- ATl/4 

Q 4 (o,t)<rc 2 



^m^iVm*) l + iV^[Mj 
1 + N^[M t T f + nfaiNm*) 

i 

1 + ^ m i a(AT m i ) + JVj^[Jl4; 



+ 



M i M k + 



< 4 r C 

i<k 







ATM* + 



AT2 



iV 2 



$(NM l ,NM k )dL lk , 



<&{NM\NM k )dL ik 
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respectively. In order to derive these inequalities we have used the assump- 
tions on the mass sizes, (3.12), and the fact that the diffusion coefficients are 
decreasing, so that they can be controlled when the masses Ml are large. 
Then, by Doob's inequality, 



pN 

(3.16) 

and similarly, 
(3.17) 



sup \A 2 (s,t)\>- 

0<s<t<T t> 
t-s<6 



<P 



N 



2 sup |A 2 (0,t)| > - 

0<i<T 5 



<^^[l + C' + C(T)], 



P 



N 



sup |A 4 (s,t)| > — 

0<s<t<T O 
t-s<8 



<4TC[1 + C + C(T)}. 



We pass to the limit N — > oo in (3.16), and conclude the proof from the 
estimates obtained in (3.13), (3.14), (3.15), (3.17) and the choice of £• □ 

4. The hydrodynamic equation. We begin with Proposition 1. 

Proof of Proposition 1. Estimate (3.1) in Lemma 1 implies that 



E Q 



sup (to A M, fit) 



<C{T) 



uniformly in M > 0. Then (2.6) follows by letting M — > oo and monotone 
convergence. 

In order to obtain (2.5), we will show that the probability measure Q 
satisfies 

r 1 uj(m) 



(4.1) 



E Q 



sup 

*(x,s)Gi 1 [' I 'x[0,r]] ^0 

ll*lll<l 



-^,fj, s ) ds 



< oo. 



Indeed, (4.1) implies that with Q-probability 1, fit(dx,dm) = vt(x,dm)dx 
for almost every t G [0, T], where v% satisfies estimate (2.5) in the statement 
of the proposition. But Q is supported on C([0,T], A4\(T x M+)), hence the 
result. 

We must therefore prove that 



E Q 



sup 

ken Jo 



to 



* ,A* S )ds 



lim lim lim E w 

K^oo A— yoo N—>oo 



sup 

l<k<K JO 



T 



oj(m) 



m 



y k > A ,Hs)ds 



< oo, 
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where we denote 

^ k ' A (s,x) = mm{^ k (s,x),A}, 

{^ fc } fceN a dense family in C([0,T],T) nB^^T] xT)] in the supremum 
norm, i?i[L 1 ([0,T] x T)] the unit ball in L 1 [[0,T] x T]. The measures /i s are 
nonnegative Q a.e., so we may take ^ k > 0, fcsN. 

For each and m > 0, let then u k ' A (s,x,m) be the solution to 



fc,A , ajm) k ,A 

Us ~r ^ U X x 

u k,A(r,-) = 0. 



fc,A 



We have the representation formula 

r-T 

.k,A 



(4.2) u k ' A (s,x,m) 



s JT 



p(a(m)(u — s),x, z)^ k ' A (u, z) dz du, 



where p(t,x,z) is the Brownian transition density on T. 
Ito's formula applied to "["^ u k ' A yields 



J o 



NM 



N 



(4.3) 



+ r T j^ u;(NM i )du k ' A d 



N 



dx 



+ f S2D N {X\NM\NM j ){u k ^)dE lj . 



if Dm now denotes 



D N {x,m,m')(f) 
1 



N 



uj(m + m')f(x, m + m ) — ui(m)f(x, m) — oj(m')f(x, m')]. 



The last term in the expansion (4.3) is nonpositive by formula (4.2), the 
assumption that > 0, and the subadditivity in m of a{m)~ l / 2 uj{m). We 
thus get 

1 ujirn) 



sup 

l<k<KJ0 



m 



(4.4) 



<E l 



EojiNm 1 ) t. a , ,• ,-. 
- Ar ; n fe ' A (x\ A/W) 



Kk<K' 



N 
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+ E* 



sup 


L 


\<k<K 





rT 



E 



uj(NM 1 ) du k 



fc,A 



N 



dx 



The second term on the right-hand side above can be easily bounded by re- 
placing the supremum by a sum over 1 < k < K and computing the quadratic 
variation of each of the resulting orthogonal martingale terms. The sum of 
these quadratic variations vanishes in the limit N — > oo; in order to see 
this, it suffices to replace u k,A by its representation (4.2) and then apply 
assumptions (2.1b), (2.3) and estimate (3.1). 

Finally, the hypothesis on , v, v* and the fact that ||\I/ fe ' A ||i < ||* fc | 
1 imply that 



< 



lim E p 

N-too 



sup 

Kk<K 



E 



u!(Nm l 



u k ' A (x\Nm l ) 



sup 

Kk<K 



m 



u k > A ,v 



< sup 

Kk<K 



uj{m) 



|^ fc,A ||i 



T(v*) <oo 



holds uniformly in A,K. 
We pass to the limit A 

□ 



oo and then K — > oo in (4.4) to obtain (4.1). 



Proof of Theorem 2. It will be enough to consider / € Cf (T x R + ) 
compactly supported and Lipschitz in m, and then use bounded conver- 
gence to obtain (2.7) for a general / S C|(T x We need to analyze the 
difference 

Z f {t) = Y,Mlf{XlNMi)-Y,m l f{x\Nm}). 



We start by writing the semimartingale Zf as 

Z f (t)=Hf(t) + A f (t), 

where T-Lf is the martingale obtained by adding the fluctuation terms arising 
from the free particle dynamics and the stochastic coagulation phenomena. 
These can be proved negligible by applying Doob's inequality, the integra- 
bility assumptions on a(m) stated in Section 1, and Lemma 1. The term Af 
is given by 

■* jaCJVArja 2 / 



dx 2 



(Xl,NMl)ds 



+ 



E 

i<3 



[M i (f(X i ,N(M i + M j )) - /(X^NM*)) 
+ M j {f(X j , N(M { + M j )) - f(X j , NM j ) )] 
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<S>(NM\NMi) Jri , 

x — — - dL 10 . 

N 

The first term of A * will clearly have the limit 
f* f /a{m)d 2 f \ , 



JTxR + 



We can guess the limit of the second term from the occupation times formula; 
we should recover the second term in the hydrodynamic equation (2.7). In 
order to obtain this expression, we replace dU 3 in the second term of Af 
by [a(m) + a(m')\V e (X l — X 3 )ds, where V e (x) approximates the Dirac 5- 
function at the origin as e — > 0. The new integral will converge weakly to 



r f(x,m + m') - f(x,m) f(y,m + m') - f(y,m') 



JTxTJM+x 



?7i m 
x k(iji, m')V t (x — y)fi s (dx, dm)/i s (dy, dm') ds 



as N — > oo. 

We will justify this exchange by showing that there exists a sequence of 
measurable sets Cjy e t with 



limlimsupP [Cn,e,t] = 0, 

€->0 jv^-oo 



such that Y7v,<E,/(i) given by 



Y^[M i (f(X i ,N(M i + M j ))-f(X i ,NM i )) 

i<j 

§{NM\NM 3 ) 



+ M 3 (f(X 3 ,N(M l + M 3 )) - f{X 3 ,NM 3 ))\- 



N 

x (dL ij - [aiNM 1 ) + a(NM j )]V e (Xi - X 3 )ds) 

satisfies 



(4.6) limlimsup.E' 



sup \r N>eJ (t)\i e c 

0<t<T 



0. 



Here C c denotes the complement of the set C. 

Suppose that (4.6) holds, and let 5 > 0, I £ N. Define 

K l (rn, rri) = «(m, m')l {l -i< m < l} l {l -i< m ,< l} , 
a\m) = a(m)l {i -i< m < /} 
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and 



sup 

0<t<T 



(Vt,f) - Oo,/> 



1 ; 1 u ,d 2 f 



J(TxR + ) 2 



2" v '"' dx 2 

f(x, m + m') — f(x, m)\ 
m' 



V e (x-y) 



x K l (m,m')fj, s (dx,drn)fi s (dy,drn')ds 



<5 



Then T^g is closed in C([0,T],Mi(T x R+)) with the Skorokod topology. 
By Proposition 1, Q almost everywhere, 



lim 

l— >oo 



7 1 u 



and 



lim 

l— >oo 



1 1 ' [f(x,m + m') — f(x, m)\ 



m 



V t (x - y)K l (m,m'),fi s <g> fj, s \ ds 



[f(x,m + m') - f{x,m)] 



m: 



V e (x - y)K(m,m'),n s fi s ) ds 



for all t€ [0, T]. These imply 



(4.7) 



lim sup lim sup Q [Fi^j] < IhnsupQI^z^a] 

/too ATfoo Zfoo 



< Q[Foc,e,2s]- 

Now, we know that fi disintegrates as fi s (dx,dm) = v s (x,dm) dx, s G 
[0,T], Q a - e - Letting e — > in (4.7), by Lebesgue's differentiation theorem, 
dominated convergence and (4.6), we have 

1 = lim lim sup lim sup Q [J-'i^s] < um Q [J r oo ) e,25] 
= 3o,2<5] 5 

if 26 is obtained replacing V^(x — y) in the definition of J-^ £j 25 by the 
Dirac function evaluated at x — y. Since 5 > is arbitrary, this implies that 
(2.7) holds with Q-probability 1. 

It remains to prove (4.6). For each e > 0, let the approximation of the 
Dirac 5 function V e be such that there exists a function u e in C 2 (T — {0}) n 
C X (T) with support contained in (—1/2,1/2), so that 

// \ 1 

IK||t»<7ie, w e(°) = 2' 
IKHt»<72, \imu' e (x) = 0, x ^0, 
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and 

u"(x) = W e (x), x^O; u"(x)<0 ifO<x<e, 

where W e is a real valued function such that | W e \ = V e , and 71 and 72 are 
positive constants independent of e. 
Consider the finite stopping time 



r e = inf|t<T,iV^[M t f >-U- 



and define 

C N ,e,T = {T e <T}. 

Lemma 1 then implies 

P N [C N ,e, T ]<C(T)V~e 

and clearly lime^ohnisupjv^^ P^Ca^t] = 0. 
Let 

G(x, y, m, m') = G N j, t (x, y, m, m') 

= [m[f(x, m + m!) — f(x, m)] 

<&(m, m') 



+ m'{f(y,m + m') - f(y,m')]]u e (\x - y\)- 



N 2 ' 

The proof follows the usual pattern after this point: we will stop the process 
upon achieving the stopping time r e , so that we may assume that P N is 
supported on 



i 

We will then apply Ito-Tanaka's formula to 



Y,G{XlXiNMlNMi) 



i<j 

in order to recover Tjv, e ,/ from the nondifferentiability of u e (|x|) at the origin 
and the particular choice of u e . One then has to check that the remaining 
terms of the expansion converge to uniformly on [0, T], when taking N — > 
00 and e — > 0, in that order. 
We have 



J2G(Xi,x{,NMi,NMi) 



i<3 



YsGix^x^Nm^Nm?) + i*r w (t) +H N , eJ (t) +A w (t), 



i<j 
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where HN,e,f(t) is a P^-martingale and Aj^^j(t) is a continuous, bounded 
variation process. We start with the former: 
ft f\r< 

i j^i 

+ E / ( E G(N[M* + M k ],NM j ) 



- G(NM\NM j ) - G(NM k ,NM j ) 



G{NM\NM k 



dE 



ik 



§(NM\NM k ) 



N 



dL 



ik 



The quadratic variation of the first term on the right above, the Brownian 
martingale, can be easily seen to vanish when N — > oo. We proceed to show 
how to bound one term in the quadratic variation Q c of the coagulation 
martingale. Consider then 

Qc(*) = /'E^K E [f(X\N{M i + M k + M^)) 



f(X i ,N(M i + M k )) 

f{X\N{M i + M j )) + f{X\ NM 1 )} 



xu e (\X-Xi\) 
<&(NM\NM k 



^(N(M i + M k ),NMi) 



N 



N 



dL lk . 



We now use that / has compact support. Let L > be such that f(x, m) = 
whenever \m\ > L. The expression between brackets in Q\ will thus vanish 
whenever NM i > L, so that we may bound Q\(T) by 

■'O ■ L 



x - 

< r (/> 2 / J2^(NM\Nm k )dL ik 
J° i<k 
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E 1 



i<k 
pN 



r 1 



< i, 



from where it follows that lim^o hmjv-s>oo -E piV [Q*] = 0- Similar considera- 
tions prove that the expectation of the rest of the terms in Q c vanish in the 
limit N —> oo, e — > 0. Doob's inequality then implies 



lim lim sup E 



(4.8) 

The Lipschitz property of / yields 



sup \n N:e j(t)\i c - Nc , 

0<t<T 



0. 



E l 



(4.9) 



]T G(Xl,xi,M l t , M>) - G{x\x\m\mi) 

i<j i<j 



The process A^^j{t) equals 

r <9 2 / 



<r(/)e[l + C(T)]. 



;2 (X\iV(M l + M^))-^(^,iVM l ) 



+ 



■ $(NM\NMi) 
x u e (\X l - ^ ^ > a(NM l ) ds 

^{X i 1 N{M i + M^))-^-(X\NM i ) 
ox ox 



x signpP - X'WcdX* -X j \) 
x — i J -a(NM l ) ds 



+ / ^[G(XSX J \iV(M i + M fc ),iVM J ') 

7o n</c jyi,fc 



G(X\X j ,NM\NMi) 
-G(X k , X j , NM k , NM j )] 

-G(X i ,X k ,NM i ,NM k ) 
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*(NM\NM k ) ik 
N 



= I 1 N,ej(t)+I 2 N,ej(t)+lN,ej(t), 

where sign(x) takes values 1 or —1 according to whether x > or x < 0. It 
is easy to see that 



(4.10) 
We can then bound 



lim lim sup E 



SU P l-Ttf.e./WI 
0<t<T 



0. 



(*)|<r(/) /^M>'j(|^-^|) 



$(NM\NMi) 



N 



a(JVAT) ds 



<r(/) 



Since lim e _>.o n^(x) = for all x ^ in T, dominated convergence implies that 



(4.11) 



lim lim sup E 



SU P l 7 Ar,e,/(*)l 
0<t<T 



as well. Finally, the fact that / is Lipschitz and (2.1a) yield 



l4,«,/WI<r(/,*)f 



^MKNMiy 
1 

iV2 



<$>(NM\NM k )dL 



ik 



so that 
(4.12) 



lim lim sup E 



sup |/L, f (t)|l 



0<t<T 



0. 



The limit (4.6) is immediate from estimates (4.8), (4.9), (4.10), (4.11) and 
(4.12), and the result follows. □ 

5. Uniqueness of the solution. In this section we seek to establish the 
uniqueness in an appropriately defined class of the solution to the hydrody- 
namic equation (2.7), which in differentiated form can be written as 



(5.1) 



1 d 2 
fit = -o(m) Vt + K(nt), 



Ho = v. 
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Here #^ denotes the partial derivative of second order with respect to the 



space variable, interpreted in the weak sense, and K is the coagulation kernel 
given by 



</,#(//)) 



[f(m + m') — f(m)] ^ — ^-/j,(dm)fi(dm') 



in' 



[(m + m')f(m + m) — mf(m) — m'f(m')] 

K(m,m') , 

x — ^(dm)fi(dm ), 

mm' 

if / is a bounded test function, /i E M f (K+) a finite measure such that n(m, 
m')/m' G L 1 (ii/x x cfyi). 

We will work under the assumption that there exists a pair of functions 
w and w bounded on each compact subset of (0,oo), such that wu^ 1 is 
bounded, u, a~ l l 2 uj and a~ l / 2 ujw are subadditive, and such that the follow- 
ing inequalities hold: 

(5.2) K,(m,m') < oj(m)u;(m'), 

(5.3) K(m,m') < uj{m)w{m') + w{m)u{m'). 

We will also require that a be nonincreasing and the initial measure v satisfy 



(5.4) v(dx,dm) <v*(dm)dx 1 
Let us define the class B(uj) by 
B(u) = \r]€ C(R+,M f (T x 



v* such that 



m 



< oo. 



sup 

t>0 



LJ(m) 



m 



„)) : rj t (dx, dm) = vt(x, dm) dx, 



< oo and sup(?n, < oo VT > }> . 



The choices 

uj(m) = [1 + c + a(l)] [m p + a(m) + 1] and w{m) = a(m) + 1 

satisfy conditions (5.2) and (5.3) in the particular situation treated in Sec- 
tions 3 and 4. In this case the coagulation kernel is given by 

K,(m, m') = <l?(m, m') [a(m) + a(m')] 

with symmetric rate verifying 



$(m, m') < c(m p + m /p )l 



[m>0,m'>0] ) 



o<p<i 
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and diffusivity a(m) such that cl~ 1 I 2 uj is subadditive. Then Theorem 2 yields 
an existence result in D(oj). The uniqueness of this solution in the larger class 
B{uj) follows from Theorem 3 below. 
Here is the main result of this section: 

Theorem 3. Assume conditions (5.2), (5.3) and (5.4) on the coagula- 
tion rate k and the initial measure v . Then for any T > 0, (5.1) has at most 
one solution {^t}o<t<T in B{u). 

We state the theorem on T to avoid introducing more terminology; in 
fact, the proof holds in M, d for a general diffusion model with coefficients 
given by a, undergoing coagulation at a rate determined by k. In that case 
we require that (5.2), (5.3) and (5.4) be satisfied by a pair of maps u,w 
such that w^uj is bounded, as before, and and a d / 2 ujw are both 

subadditive. 

The result will follow from considering an approximating system of equa- 
tions to the coagulation equation, for which existence and uniqueness can 
be easily derived. The method is a simplification of a technique developed 
by Norris in [16], we have thus tried to adhere to his notation whenever 
possible. We are able to make a significant shortcut in the proof due to 
the assumption that we already have got one solution in B(u), this yields 
a crucial monotonicity property in the approximating scheme as a direct 
byproduct of its construction (compare Lemma 5.1 below with Lemmas 5.5 
and 5.6 in [16]). 

Before we can proceed to prove the theorem we need to introduce some 
definitions. Given s > and /x G Mf(T x R+), we will denote P s fi G Mf(T x 
M+) the measure given by 

(f,P s fi} = J f(z,m)p(a(m)s,x,z)fi(dx,dm)dz, 

where, as before, p(t, x, z) is the Brownian transition density on T. Hereafter, 
we will use the abridged notation p^' z (m) to denote p(a(m)t,x, z). We also 
introduce the kernels K + and K~ on ,M/(]R+) defined as 

K+^)(dm) = [ K(m ' '™"\ (dm')ii{dm'% 

J m'+m"=m ™ 

K-(n)(dm) = [ Kijn,m '\ {dm)n{dm'). 
Jr+ rnf 

With this notation, we now show that (5.1) is equivalent for /i £ B(u) to the 
integral equation 

(5.5) Mt = ^+ / Pt- r K + (fi r ) - P t . r K-(fi r )dr. 

Jo 
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By integrating against a test function and differentiating in time, it follows 
that any solution to (5.5) satisfies (5.1). Conversely, let \it be a solution to 
(5.1), and set fit equal to the right-hand side of (5.5). Then n — jit verifies 

d Id 2 

— (jlt-fit) = -a(m)-^(jlt-Ht), /2 -Mo = 0, 

and we conclude that jit = lit- In particular n t is a solution to (5.5). 

Given a bounded map c s (x) : [0, T] x T — > K, let p m (c) be the propagator 
associated with the operator ^o(m)^^ — c s (-) on T, and consider the kernel 

Pts(c) n(x , dm) = J /j,(dz,dm)p m (c)(s, z;t,x) dz, 
Jt 

fi in Mf(T x R + ). 

Let E n = [l/n,n] for n£i, and define , if~ by analogy with ET + : 
K+(/j)(dm) = K + (n)(dm)l {meEn} , 

K-{n){dm)= / - t h m+m 'eE n \^{drn)^{dm') 

JR+ m 

and K n = K£ — K~ . Define v n as the restriction of v to E n , u n (dm) = 
v{dm) x l{ meEn }- 

In order to keep the number of definitions to a minimum, given a test 
function g, we will use {{g,^)) to denote the integral {{g,v)) in the case 
that the decomposition fi(dx, dm) = u{x, dm) dx holds, where we recall that 
{(g,v)) stands for the single integral J g{m)v(x,dm). 

The proof of Theorem 3 will follow from the following lemma. 

Lemma 3. Let {h s }o<s<t & £>( w ) be a solution to (5.1) with initial value 
v, and assume that conditions (5.2), (5.3) and (5.4) hold. Then, for each 
n£N, there exists a unique kernel (jlt)o<t<T m S(uj), such that 

(5.6) A? = ^(c n K+ f P ts {c n )[Kt(K) + ^K]ds, 

Jo 

where 



cUx,m) = (L(m)^,P t u\\ + J o (L(m)^,P t - s [K n (jl-)}y\ds. 

Moreover, jl™ satisfies 

K<K +l <^s for all s€[0,T\. 
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Proof. The method of the proof is classical. We will take advantage of 
the fact that all relevant quantities are bounded in E n to define a Picard 
iteration procedure, and then prove by a contraction mapping argument that 
the scheme converges to a solution, which is finally shown to be unique. 

Define 



-t,o 



^ m )^,^)) + Jf'^ m) ^,p ( _, [ i r „ K , )])) (is . 

Condition (5.4) on v and the fact that p, n is supported on E n imply that 
c™ (x,m) is well defined and bounded. 
'Let 

c t (x,m) = (Au!(m) ,nt ) )■ 

We claim that < ct(x,m) < c^q" (x,m) < cf (x,m) < Cuj(m), where C is 
a positive constant that depends solely on /i. Note that pu' x (m)u>(m) is m- 
subadditive, which can be easily seen by expanding pu' x as a series and using 
the subadditivity of and the monotonicity of a. We next write 



oj(m' 



\Pt- s ( m ' + Tn")u)(m + m") 



m 

1 

y.x i i\ ( i\ y,x ( ii\ ( 

-PtL s {m )u(m ) -p y t L s (m )u(m )J 

x K{m ';™P lE n (m' + m'')ft )0 (y, dm'^y, dm") dy. 
mm 

Together with the subadditivity of pu' x (m')to(m'), the inequality j2™ < 
ft™ o 1 — ^ s anc ^ the obvious inclusion E n C E n+ \ , this implies 

(5.7) 

> c™ (x, m) > c"q 1 (x, m) > ct(x, m) > 0, 

as required. 

We will need the fact that the solution ji. satisfies 

Ht = P t (c)u+ [ P ts {c)[K + {fi s ) + 5 s ^ s ]ds, 

(5.8) J ° 

1 1 ijj(m)uj(m') — K(m,m') 
d t {x,m) = { { ,n t 
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This can be proved by an argument similar to the one used to show the 
equivalence of (5.5) and (5.1). 

Fix fcsN, and assume that a sequence of measure paths and mappings 
fis,k - Z^fe 1 - P* an d c™ k (x,m) > c££ (x,m) > c s (x,m) > have already 
been defined. Let 

1 uj(m)uj{m') — K(m,m') 



Of 



m' 



~n 

>Pt,k 



so that < 5™ k < 5™ k+1 < 5 t . Finally, set 

Kk+i = Pt{c n k )v n + f P ts {c n k )[Kt{Kk) + ds, 
Jo 

cl k+1 (x,m) = ^(m)^,P t ^ 



The Feynman-Kac formula (cf. [17], Chapter 8) and identity (5.8) then 
imply 

0<# ifc+1 <£$!i^« 
and thus, by the same arguments that yield (5.7), 



( m )/ ^ m ' ) v*\>c n >c n+1 

\ m ) \ —/ ' V J ^ C s,k+1 C s,k+1 



> c, >0. 



The kernel p n sk is supported on E n for each n and all k, s. Also, p™ k < /j, s 
and the fact that uj{m)/m is bounded below in E n imply that the marginal 
density dp™ Adx ,M+) / dx G L°°(dx) and its L°°-norm can be bounded uni- 
formly in k, 

dp™ k {dx,R + 



(5.9) 



dx 



for all k G N. 



In view of these observations, we define a norm in the vector space 
r] G M f (T x E n ), ^(dx,R + ) G L°°(dx) 
by 

lllp- q\\\ = WWpx - &IIIL, 

if p(dx,dm) = p x (dm)dx and a similar formula holds for g. Here ||p x — 
is the total variation norm of the signed measure p x (dm) — g x (dm), which 
we may compute as 



\Px 



\p x ~ Qx\(dm), 
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if \p x — Qx\(dm) denotes the total variation of p x (dm) — g x (dm). Equivalently, 
g\\\ =sup\ (f,p- g)\f: / \ supf(x,m)\dx <l\. 

L JJ m ) 

Let f(x, m) be such that J T | sup m f(x, m)\ dx < 1. We obtain 

= (f,P t {ct)v n -P t {ct_ 1 )v n ) 
(5.10) t 

+ f (/,A s (cj?)[^+(/i? )fc )]-A s (cLi)[^(^ fe -i)]) ds 

J 

Jo 

We apply the Feynman-Kac formula to bound the first term on the right- 
hand side of (5.10) by 



TxR 4 



\clk(Xs,m) -Cs ik -i(Xs,m)\ds )f(xt,m) 



o 



v n (dx, dm), 



where Xs is a Brownian motion in T with diffusivity a(m). We have 



\ c s,k(Xs,m) - c^^ixs^ds 



o Jo 



oj(m' + m" 



ml + to 
k(to', to") 



^rfj*- (to' + to") - K) 



TO 



l{m'+m"G_E Tl } 



- AS,fc-i ( s '> dm ')Vu,k-i( x> > dm ")\ dx' du ds. 

Note that in all these expressions the mass variable takes values to, to', to" 
or to' + to" belonging to -E n , we may therefore replace all functions that 
have it as an argument by an upper or lower bound, as necessary, that do 
not depend on the mass variable. 

Now, if p and g are two finite, positive measures on (X, Q), then p® p 
and g® g are finite, positive measures on (X x X, Q x f2) and 

\p®p- g<8> Q\ < [p{X) + g(X)]\p - g\. 

In particular, 

\i^ k (x\dm')pZ ik (x',dm'') - j% k _ 1 (x\dm')i% )k _ 1 (x' 1 dm'')\ 
<2(3 n \pl k (x',dm')-pl k _ 1 (x',dm')\, 
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fi n as in (5.9). 
We thus have 



(/^(cjJK-i^K) 

<T / / p(s — u,x' ,w)p(s,x,w)p(t — s,w,y) 
Jo Jo 

x ( sup|/(y,m)| ) \\KJx') - Km-i( x> ) 



x v n (dx,dm) dx' dw dy duds, 

where T is a positive constant that depends on n and fi n and can therefore 
be chosen uniformly in k. 

Replace now v n by its upper bound is*(dm)dx, take the L°° norm of the 
total variation factor lA™^^') ~ iff k—i( x ')\i an( ^ integrate in x,x',w and s, 
to obtain 

(/, P t (c£K - Pt(c n k ^)v n ) < r(T) f \\\K, k - K >k -i\\\ du. 

Jo 

In this last step we used that J" T |sup m f(x, m)\ dx < 1. 
Similar computations yield 

Jo 

<r(T) fWu,k-K,k-i\\du 

Jo 

and therefore, by (5.10), 

</,&Vi - Kk) < T (T) f \\\K,k ~ K,k-l\\ du, 

Jo 

r(T) uniform in k, f. Take the supremum over / to conclude that for t <T, 

du. 



(5.ii) \\\ft M1 - ft tk \\\ < r(T) / p™ fc - I 

JO 



Hence, by a standard contraction mapping argument, /}J n converges in 
Mf(Tx E n ), uniformly in t < T. The limit is a continuous map fi n : [0, T] — > 
Mf(T x E n ) that satisfies (5.6). Moreover, the properties that \i n sk < fi^ 1 < 

fi s and c™ fc > c™"^ 1 > c s for all k,n G N 2 imply that the same holds for /2 n 
and c n : 

A"<A" +1 <At s and 4 l >c" +1 >c s , 
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for all n £ N. 

Suppose now that jl n and fj n are two solutions in V(K) to (5.6) with 
respective values of the mapping c defined as in the statement of the lemma. 
A careful revision of the arguments leading to (5.11) shows that the proof 
goes through verbatim if we replace /ijfc+i and ^ by fixed points of the 
iteration scheme fi n and fj n . We thus have 

|||A?-^lll<r(r) fWu-r u \\du, t<T, 

Jo 

which implies /i n = ff 1 by GronwalPs lemma. □ 

Proof of Theorem 3. Set Aq = lgci/, and define the kernel K~((i) = 
K~(n) — K~(n). Given mappings jl n and c n as in Lemma 3, define 

= P t v n + f P t -.[K+(jfi) - K-(j%) - 7 M ds, 
Jo 

\? = P t \%+ f P t -s[K-{K) + l n sK]ds, 
Jo 

where 

7™(x,m) = c^(x,m) - (( u{m) U , ) ) > c™(x,m) - c s (x,m) > 0. 



The inequality in the last line follows from /2™ < fi s and the definition of c"; 
it has the consequence that A" is a positive measure for all s <T,n 

We claim that fi n = ft 1 . Indeed, differentiating in the equations satisfied 
by ji and jl shows that both maps verify the equation 



2 



m = \a(m) ^ Vt + K+(ft) - [c? - W 

with initial value v n , so that their difference is a weak solution to r) = 
^a(m) -7^2 rjt started from the zero measure, and therefore it must be the 
null measure. In particular this implies that fj,™ < fi s . 
We have, by the definition of c™, 

c?(x,m) = (L(m)^,P t u\\ + jf (L( m )^,P t _ s [K n (fi-)}\\ ds 
m / / \ \ m' 



and on the other hand, from the definition of 7, 



f 
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Hence 7 »* = ((a;(m)^,A?». 

From this point on, the proof is copied from that of Theorem 5.4 in [16]. 

Let us define af = , A")). Then, due to the subadditivity of w(m)a(?n)~ 1 / 2 , 

we get 

for all n £ N, and it follows from Lemma 3 that a™ > We can hence 

define the monotone limits 

a t = lim a", fi t = lim fj,?, 

n— >oo — n— >oo 

which satisfy 

/ / u(m) \\ 1 1 (jj(m) 

Since ui(m)/m > 0, in the case that at = a.e. we can conclude that \it = 
a.e., for all t £ [0, T]. The uniqueness of the solution to (5.1) will then follow 
from the uniqueness of the solution to (5.6). 

We will now show that at vanishes. We start by proving that h n (t) = 

sup s>0 ||(( a; , PsfJ-t)) Woo can be bounded uniformly in n and t E [0, T]. We 
apply P s to the definition of fif, multiply by w 2 (m)/m and integrate over 
E n to obtain 



m II \ \ m 

(5.12) 

By the subadditivity of lj and wa^ 1 / 2 , for any u > 0, x, z 6 T, we have 
u) 2 (m + m')p^ x (m + m!) - u 2 (m)pl' x (m) - oj 2 (m')p z ^ x {m) 
< u(m)pu' x (m)u}(m') + uj(m)u;(m')pl' x \m') . 



Then 



l \p s+ t^KM)))(x) 



in 



<2 ( u{m)p z s ,x t _ r {m)oj{m')— — - - pL^{z,drn)^{z,drn')dz 



< 2 
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UJW 



m / / \ \ m 



00 



< 2 



U1W 



in 



u n 



tu- 



rn 



Ps+t-r^r ) ) 0) 



+ 2 



_ r 

m 



— ,P s+ t- r tf)){x), 

m 



where we used the bound n(m,m') < uj(m)w(m') + w{m)ui{m'). Now, if we 
replace ui 2 by uivo in (5.12), the subadditivity of a~ l / 2 uiw implies that the 
time integral term is nonpositive, and therefore 



sup sup 

r s>0 



UJW 



m 



< ( , V 

m 



< oo 



by condition (5.4) on v. This assumption also implies that 



sup 

s>0 



m 



<r<oo 



for some positive constant V. Taking the supremum over s > in (5.12), we 
conclude that 

K(t) < T + 2T' [ h n {r) dr with V > 0. 
Jo 

Note that the constants T, V can be chosen independently of re. Then h n (t) < 
Te 2T T holds uniformly in n and t <T, as claimed. 

We now consider the L 1 norm of af. We replace 7™(x,m) by its upper 
bound uj(m)at(x) in the definition of A n and pass to the limit as n — > oo. 
By dominated convergence we have 

ui 2 (m) 



Mil < 



< 



JTxTxR + 

sup sup 

n s>0 



m 



-a s (z)p t , _ s (m)fi (dm,dz) dx ds 



o 



CO 



m 



\ a s\\i ds 



< r(r) f 

Jo 



\oi s \\\ ds. 



Since ||at||i < {u)/m,v*) < oo, this implies at = a.e., as required. The the- 
orem follows. □ 
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